1 



NNLO unquenched calculation of the b quark mass. 
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By combining the first unquenched lattice computation of the B-meson binding energy with the recently cal- 
culated two-loop contribution to the lattice HQET mass, we determine the MS b-quark mass, m(,(mi,) at the 
NNLO. We find mt(mb) = (4.26 ± 0.03 ± 0.04 ± 0.05) GeV. The inclusion of the two-loop effects is one of the 
steps necessary to extract rnt(rnt) with a precision better than 0{Aqcd), which is the uncertainty due to the 
presence of an IR renormalon singularity in the perturbative series of the residual mass. Our results have been 
obtained on a sample of 60 lattices of size 24^ X 40 at /3 = 5.6, using the unquenched Wilson action with two 
degenerate sea quarks. The quark propagators have been computed using the unquenched links generated by the 
TxL Collaboration. 



1. Introduction and Motivation. 

Quark masses are fundamental parameters of 
QCD which cannot be determined by theoretical 
considerations only. Moreover, they are very im- 
portant for phenomenology because enter many 
theoretical predictions of physical quantities like 
the CKM matrix elements, the CP violation pa- 
rameter e'/e, the AJ = 1/2 K ^ tttt amplitude, 
the B meson semileptonic decays, the B — B mix- 
ing . . . etc. Quark masses cannot be measured 
directly since quarks are confined in the hadrons. 
Therefore, a short distance definition of the quark 
mass, which is scale and scheme dependent, must 
be adopted. 

In this paper, we report on the first unquenched 
HQET lattice calculation of the b quark mass. 
The main idea is to combine unquenched 
HQET lattice computations of the B-meson prop- 
agator with recent NNLO analytical formulae of 
the matching of the continuum A'lS quark mass 
to the lattice one ||]. We stress that both the 
unquenched lattice simulation and the NNLO 
matching are necessary ingredients to improve the 
accuracy of the results: the former is necessary 
to control potentially large vacuum polarization 
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contributions to the B meson propagator and the 
latter is crucial to cancel the renormalon ambigu- 
ities in the matching 

2. How to determine the b-quark mass on 
the lattice. 

The key idea is to match the 6-quark propa- 
gator in QCD to its lattice HQET counterpart. 
Since the lattice HQET is an effective theory, the 
relation between these propagators is, to lowest 
order in 1/mt,, 

S^^{p,mb;^i) = C{mb/n,as) Sj^'^dv ■ k)a;^i) (1) 

where p = rrib v + k is the external momentum of 
the 6-quark, v is its velocity, k is the residual mo- 
mentum with ^ mb, fJ. is the renormalization 
point, a is the lattice spacing and C{mb/ ^,as) 
is the Wilson coefficient. Note that, as expan- 
sion parameter of the HQET, we have chosen the 
quark mass m^. In this way, all the mass depen- 
dence is factorized in the coefficient function C. 
The method to fix C is well known: it consists in 
calculating the 6-quark propagator in QCD and 
in the lattice HQET to a given order in as, com- 
paring both expressions at a fixed jj, = /jq and ex- 
tracting C{mb/fJ.Q). Renormalization group can 
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then be used to evolve this function to lower 
scales. The important point is that by rewrit- 
ing eq. (|l]) in terms of the pole mass, nipoie, is 
easy to find the relation 



Y 

T-pole = TOfc + > {as{a)) 

^-^ a 



(2) 



where the last term is the perturbative expan- 
sion of the residual mass generated in the lattice 
HQET, Sm. The coefficients Xn arc functions of 
\n{mba). 

On the other hand, the HQET mass formula, 
to lowest in I /mi,, 
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£ ^ rrib + 0{l/mb) 



(3) 



allows us to write the unknown mass mt in eq. (g) 
in terms of the physical B-meson mass, Mb- The 
so-called binding energy, £, is independent of rn^ 



mpoie = Me 



£ 



n=0 



(4) 



£ is not a physical quantity because it diverges 
linearly as a — > 0. Since the pole mass is finite in 
this limit, this divergence is cancelled by the last 
term of eq. . In practice this cancellation is not 
perfect and hence one cannot take a too small. 
As a large nj calculation demostrates (see [||), 
the series in eq. (||) has IR renormalon singulari- 
ties, the same as the pole mass. In other words, 
the coefficients Xn grow as n! as n ^ oo. This 
behaviour gives rise to ambiguities of 0{Kqcd)- 
In order to avoid this problem, it is convenient 
to use a short distance definition of the b-quark 
mass, like the MS, mb, which is free of IR renor- 
malons. The relation between the pole and the 
MS mass is known to 0{a'l) [Q and has an IR 
renormalon that cancels the one in (|^) . The can- 
cellation, however, is delicate so that one need 
to know many coefficients X„ to keep under con- 
trol the partially removed renormalon ambiguity. 
Xq is obtained easily in terms of a one-loop lat- 
tice integral. Recently, Martinelli and Sachrajda 
have performed the calculation of the coefficient 
Xi |0] . For the Wilson action, they are 



Xo = 2.1173 

Xi = (3.707 - 0.255 n/) ln(TO6a) 
- (1.306 + 0.104n/) 



(5) 



Putting all together, we find 

/ 2.1173 

mbljflb) = Mb - £ ^ asirrib) 



+ - [3.197 In(TOfca) - 1.514] asimb)' 



4 as{mb) 



9.58 



asirrib) 



(6) 



where we have taken n/ = 2, the number of sea 
quarks in our simulation. Note that the remain- 
ing ambiguity in our result is only ©(Aqp^/to;,) 
which is beyond our precision. By computing 
£ through lattice simulations of the HQET and 
by using the experimental values of the B-meson 
masses, we can extract nib. 



3. Lattice computation of £. 

As well known, £ can be extracted by studying 
the large time behaviour of the two-point func- 
tion of the B meson, C2 in the HQET, 



C2(t) 



z 



-£t 



(7) 



The light quarks are described by the Wilson 
action. In order to improve the isolation of 
the ground state, we use cube and double cube 
smeared axial currents as interpolating operators 
of the B meson Q] . The actual value of £ has some 
dependence on the cube size and the smearing 
type due to contamination from excited states. 
To obtain our best estimate, we compare differ- 
ent methods and account this systematic effect in 
the final error. We use the Standard Method, in 
which we base our results on the best cube, de- 
fined as the one which yields the largest and flat- 
test effective mass plateau. In practice, we have 
few cube sizes (only two in our simulation, 7 and 
9), and hence the best cube is difficult to find. 
To improve our results, we also use the Multifit 
Method which consists in performing a global fit 
of the data for all smearing types and cubes sizes 
imposing that the binding energy be the same for 
all of them. In order to reduce the effect of ex- 
cited states, it is convenient to fit the data to a 
two-state form of C2 rather than to eq. (Q) . 
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4. Analysis and results. 

We performed an unquenched Wilson simula- 
tion with two degenerate sea quarks at two values 
of their mass, fc^ea — 0.1575 and kgea = 0.1580, 
at /3 = 5.6 on a 24^ X 40 lattice. The sample 
is of 60 configurations at four values of the va- 
lence light quark masses. We compute the quark 
propagators and correlation functions from the 
gluon configurations generated by the TxL Col- 
laboration. Since our simulation is unquenched, 
the procedure to measure the lattice quantities, 
is slightly different from the usual quenched case. 
Since there is some confusion in the literature, we 
want to stress a point that we consider important 
in unquenched analyses. We think that the cor- 
rect strategy consists in performing an indepen- 
dent quenched-like study of all lattice quantities, 
including the lattice spacing, for each ksea- Only 
at this point, the quantities expressed in physical 
units can be extrapolated as a function of the sea 
quark masses. The reason is that a change of the 
sea quark mass modifies the value of the coupling 
constant, and hence, may induce a rapid varia- 
tion of the value of the lattice spacing. There- 
fore, lattice results from different values of kgea 
are not directly comparable until they have been 
converted to physical units. We find, from mK* , 

«"'Uk. = {2.51(6), 2.54(6)} GeV (8) 

and for the binding energy 

aSs^ ^ {0.588(11)(5),0.606(15)(2)} 

a£B, = {0.620(8)(4),0.632(12)(2)} (9) 

for the two values of kgea- More details will be 
given in 

We have carefully studied different sources of 
systematic errors in eq. (|^). We find that the 
dependence of mb(mb) on the unknown Ag^j^ 

is small. Since the quenched value is Ag^j^ ^ 
250 MeV 1^, and the physical one is expected 
to be larger, we varied the value of Ag^^ in the 
range [250, 350] MeV. We also tried other options 
to evaluate a"^"^ In the numerical evalua- 
tion, we also used eq. (^) expanded to 0{al). In 
this way, we obtain a rough estimate of higher- 
order terms which is found less than 3%. Chi- 
ral extrapolations are also under control because 



the b~quark masses obtained from the Bd and 
Bs mesons are nicely compatible. Results from 
different smearing methods are also compatible, 
showing that the ground state has been well iso- 
lated. Finally, the dependence on fc^ea is very 
small. This allows us to take for our best esti- 
mate of rribijnb) the value obtained at the light- 
est one, ksea = 0.1580. Taking into account all 
errors, we find 

mf,(TO6) = (4.26 ±0.03 ±0.04 ±0.05) GeV (10) 

where the first error is statistical, the second is 
systematic and the third is an estimate of the ef- 
fects of higher-order corrections in the matching. 
It is interesting to compare our result with our 
old quenched numbers at three different lat- 
tice spacings = 2.0(2), 2.9(3) and 3.8(3) GeV, 
which have been reanalyzed in |^ to include the 
two-loop correction to the matching 

mf,(TOfc) = {4.34(5), 4.29(7), 4.25(7)} GeV (11) 

where the error comes from the lattice uncertain- 
ties. In this case, an educated estimate of higher- 
order corrections is 0.10 GeV. As can be seen, 
they are in good agreement with our unquenched 
result. Given the present uncertainties, in partic- 
ular the small physical volume used to obtain the 
result for — 3.8 GeV, we are not in the po- 
sition to attempt the extrapolation of the results 
in (Inl). 
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